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Abstract 

Observations from the field and experimental evidence suggest that different strains of infectious salmon anaemia 
virus (ISAV) can induce disease of varying severity in Atlantic salmon. Variation in host mortality and dissemination 
of ISAV isolates with high and low virulence was investigated using immersion challenge; from which mortality, 
pathological, immunohistochemical and preliminary molecular results have been previously published. Here, 
real-time RT-PCR analysis and statistical modelling have been used to further investigate variation in virus load and 
the response of four select immune genes. Expression of type I and II interferon (IFN), Mx and ylFN induced protein 
(ylP) to high and low pathogenic virus infection were examined in gill, heart and anterior kidney. In addition, a 
novel RNA species-specific assay targeting individual RNA types was used to investigate the separate viral processes 
of transcription and replication. Unexpectedly, the low virulent ISAV (LVI) replicated and transcribed more rapidly in 
the gills compared to the highly virulent virus (HVI). Subsequently LVI was able to disseminate to the internal organs 
more quickly and induced a more rapid systemic immune response in the host that may have offered some protection. 
Contrary to this, HVI initially progressed more slowly in the gills resulting in a slower generalised infection. However HVI 
ultimately reached a higher viral load and induced a greater mortality. 



Introduction 

Infectious salmon anaemia virus (ISAV) is an aquatic or- 
thomyxovirus that infects Atlantic salmon. The virus 
can cause serious disease problems resulting in substan- 
tial financial implications for the aquaculture industry. 
Infectious salmon anaemia (ISA) has occurred annually 
in Norway [1] since its discovery in 1984 [2] and has 
also been reported in North America [3-5], Scotland 
[6,7], Faroe Islands [8,9] and Chile [10]. ISAV is envel- 
oped and possesses a genome of 8 negative-sense single- 
stranded RNA segments, coding for at least 10 proteins, 
similar to that of influenza virus [11,12]. Segments 1-4 
produce the proteins that comprise the viral ribonucleo- 
protein core. Segments 5 and 6 encode the fusion (F) 
[13] and haemagglutinin-esterase (HE) [14,15] surface 
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glycoproteins respectively. Segment 7 possesses two open 
reading frames (ORFs), the first of which is thought to 
produce a minor or non-structural protein (NS) [16] that 
interferes with interferon (IFN) and IFN induced systems 
[17,18]. The second ORF, generated via splicing, is be- 
lieved to encode the nuclear export protein (NEP) [19]. 
Segment 8 also contains two ORFs; ORF1 encodes the 
major structural matrix (M) protein [14,16] and ORF2 is 
believed to produce a minor structural protein also with 
IFN antagonising activity [17]. The M protein is the most 
abundant protein in the virion [20] and therefore RNA of 
this segment is readily targeted in diagnostic tests [21]. 

The pathogenicity of ISAV is a multifactorial trait 
dependent on the function of viral proteins, interactions 
with host immune responses as well as various environ- 
mental factors. Little is known regarding the specific 
interplay of these aspects and how they influence ISAV 
infection. Nevertheless there are multiple reports that 
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different strains of ISAV display different levels of viru- 
lence and associated mortality [22-26]. A deletion within 
the highly polymorphic region (HPR) of segment 6 has 
been suggested to be an important virulence marker 
[25,27]. All pathogenic strains reported to date contain 
an HPR deletion with respect to the putatively ancestral, 
non-virulent HPRO virus [9,28]. In addition, substitu- 
tions and insertions adjacent to the putative cleavage site 
in the F protein have also been linked to virulence 
[27,29]. However, even isolates with identical HPR and F 
protein sequences have shown variation in virulence 
[25,27] indicating additional factors are also important. 

Vertebrate nucleated cells rapidly secrete the cytokine 
molecules interferon (IFN) in response to viral infection. 
In turn, IFNs activate complex signalling pathways indu- 
cing the expression of hundreds of direct and indirect 
antiviral genes [for review see [30,31]]. In fish, both the 
Type I and Type II IFN systems are now known to be 
imperative in the innate and adaptive antiviral responses 
respectively [31,32]. ISAV is a potent inducer of immune 
genes yet the response offers little protection [25,33,34]. 
Intraperitoneal injection of a highly virulent ISAV indi- 
cated significant monophasic induction of several im- 
mune genes in kidney on day 6 post-infection, followed 
by increased virus production and high mortality [33]. 
The antiviral protein Mx has become a direct indicator 
of the innate Type I IFN response, induced by IFNs a 
and (3 in a variety of cell types. Similarly, the IFNy- 
induced protein (yIP) gene can be used as an indicator 
of the adaptive Type II IFN response, stimulated by ylFN 
specifically in immune cells. The expression of these 
interferon stimulated genes (ISGs) as induced products 
of each system showed both the innate and adaptive re- 
sponses were stimulated concurrently by ISAV injection 
[33]. More recently, cell cultures infected with ISAV of 
high or low virulence have indicated variations in im- 
mune response between strains [35,36]. 

The present study utilises an immersion infection of 
Atlantic salmon with two ISAV strains previously classed 
as low (LVI) and high virulence (HVI) [23,37]. The mor- 
tality, pathology, immunohistochemistry (IHC) and ini- 
tial virus segment 8 qPCR results observed following the 
challenge are presented in McBeath et al. [37]. Here, the 
molecular results are analysed in more detail using stat- 
istical modelling of the qPCR data of viral segments 7 
and 8 obtained from gill, heart and anterior kidney sam- 
ples. The separate viral processes of replication and tran- 
scription were also investigated individually for both 
viruses using novel RNA species-specific assays based on 
a recently published method [38]. In addition, the ex- 
pression levels of four host immune markers were also 
monitored and statistically analysed to determine if in- 
fection with LVI or HVI had variable effects on the im- 
mune system. Ultimately, LVI was shown to replicate 



more rapidly in the gills resulting in more rapid dissem- 
ination and immune response throughout the host, yet 
HVI reached higher viral loads, causing a more serious 
infection and inducing greater mortality. The results of 
this study contribute to our understanding of the patho- 
genesis of this important pathogen in Atlantic salmon 
aquaculture. 

Materials and methods 

Virus propagation, immersion challenge and organ 
sampling 

Details regarding the virus production, infection challenge 
and organ sampling procedure have been described previ- 
ously [37] . Briefly, fish were challenged by immersion with 
either 10 4 TCID 50 of highly virulent ISAV Glesvaer 2/90 
(HVI) or low virulent Can/F679/99 (LVI), or mock- 
infected with virus-free cell culture medium as uninfected 
negative controls. After 2 h, fish were transferred to 
sampling tanks (n = 75) or to observation tanks where 
cumulative mortalities for each treatment were observed 
(n = 20). Four fish were sampled 6 h post infection (pi) 
and on days 1 to 8, 10, 12, 14, 19 and 23 pi from each 
challenge group. Organ samples from gill, heart, and 
anterior-kidney, were collected in 1 mL RNAlater (Qiagen) 
and stored at -80 °C. 

Viral and immune gene expression real-time RT-PCR 

RNA was extracted from 5 mg organ samples (gill, heart, 
anterior kidney) using the QIAsymphony" RNA robotic 
system (Qiagen) according to the manufacturers' proto- 
col and eluted in 100 uL RNase-free dH z O. The RNA 
was reverse transcribed to cDNA using the TaqMan" Re- 
verse Transcription Reagents kit (Life Technologies, UK) 
with oligo-d(T)i 6 as described previously [33] using a 
20 uL total reaction volume. Real-time RT-PCR assays 
were performed as described in McBeath et al. [37]. Pri- 
mer and probe sequences (Additional file 1) for the as- 
says targeting viral gene (no differentiation of ORFs in 
either assay) segments 7 (seg7) and 8 (seg8) from both 
European and North-American genogroups, immune 
markers Type I (al and a2), Type II interferon (IFN), 
Mx, ylFN-induced protein (yIP) and endogenous control 
elongation factor la (ELF), have been described previ- 
ously [21,33]. Absolute quantitation of transcripts was 
carried out. The cycle crossing point (Cp) values were 
converted into expression values normalised against the 
reference gene, ELF, using the statistical standard curve 
method to produce relative expression ratios [39]. 

Statistical analysis 

The analysis focussed on modelling the scale of virus loads 
and immune responses for each ISAV isolate throughout 
the experiment. Statistical analyses were performed within 
the R statistical environment version 2.15.2 [40]. 



McBeath et ah Veterinary Research 2014, 45:83 
http://www.veterinaryresearch.Org/content/45/1/83 



Page 3 of 1 3 



All observations greater than zero were transformed to 
logio, converting the analysis into an evaluation of the 
variation in the scale of detectable expression for defined 
genes associated with two ISAV strains over time. The 
distributions of the transformed observations for each 
gene-strain-time combination are also more satisfactory 
approximations of the normal distribution. A prelimin- 
ary inspection of the transformed observations suggested 
that there were differences in the dispersion of the 
scaled values at different times, and weights comprising 
the reciprocal of the variance for each strain-time com- 
bination were therefore calculated. These weights were 
used for the modelling of the changes in the scale of de- 
tectable gene expression for each strain over time using 
a locally- weighted running-line smoother [41]. Predicted 
values with their 95% confidence intervals (CI) were 
plotted. Plots exhibiting minimal overlap of CI for the 
two strains indicate differences in the trajectory of de- 
tectable expression for that gene over time between the 
two strains. A more formal evaluation of differences be- 
tween strains was carried out by estimating the reduc- 
tion in the residual sum of squares of nested linear 
models [42] describing gene expression on the potential 
explanatory variables of strain and time; /"-values of < 0.05 
were categorised as statistically significant. 

The association between the difference in the level of 
expression of the immunological markers (e.g. Mx) of 
the two strains with the difference in the load of the two 
strains (evaluated using both the segment 7 and segment 
8) was evaluated for each tissue using Spearman's rank 
correlation coefficient (r s ) [43]. The statistical signifi- 
cance of associations is indicated by "star value" (*0.01 < 
p < 0.05, **0.001 <p< 0.01, ***p < 0.001). 

RNA species specific real-time RT-PCR 

RNA species specific reverse transcription (RT) was per- 
formed as described previously [38] using tagged primers 
for mRNA and cRNA (Additional file 1). Reverse tran- 
scription was performed at 65 °C using the Thermoscript 
RT kit (Life Technologies) according to manufacturers' 
protocols. 

Real-time PCR was carried out using SensiFAST Sybr 
no-rox Master Mix (Bioline) and primer pairs consisting 
of one primer specific to the tag portion of the RT- 
tagged primers and the other specific to the appropriate 
viral sequence (Additional file 1). Reactions consisted of 
1 x SensiFAST Sybr no-rox Master Mix, 10 uM each pri- 
mer, 1 uL cDNA template in a 20 uL volume and were 
subjected to a 95 °C for 2 min pre-incubation, 45 cycles 
of 95 °C for 5 s, 62 °C for 10 s, 72 °C for 10 s, followed 
by melting curve analysis (95 °C for 5 s, 65 °C for 1 min, 
slow heating to 97 °C and cooling at 40 "C for 30 s) 
using a Roche Lightcycler LC480. The cycle crossing 
point (Cp) values were cross checked with melt curve 



data and samples were deemed positive if the melt curve 
temperature was within ± 0.5 °C of the expected value. 
Samples with no Cp value or a Cp value correlated with 
an incorrect melt curve temperature (i.e. indication of 
non-specificity) were classed as negative. Samples with 
dual melt curves were classed as positive if one was at 
the correct temperature. Due to the complication of 
non-specificity, no attempt was made to quantify the dif- 
ferent RNA types. 

Results 

Experimental infection and mortality 

In fish infected with HVI, mortality began on day 13 pi 
and the cumulative mortality reached 100% by day 23 pi. 
In the LVI group, mortality began on day 17 pi and cu- 
mulative mortality reached 20% by day 23 pi. All dead 
fish tested positive in kidney for ISAV segment 8 by real- 
time PCR. No mortality was observed in the uninfected 
negative control fish. Further details of the above results 
including pathology and immunohistochemistry (IHC) 
are presented in McBeath et al. [37]. 

Viral kinetics profile by detection of ISAV segment 8 

All infected fish tested positive for segment 8 in the gills 
at every time point. In heart, no virus was detected at 
6 h pi, 2/4 fish tested positive for LVI and 1/4 for HVI 
on day 1 pi. After 2 days, heart was positive in 3/4 fish 
infected with LVI and 2/4 infected with HVI. Kidney 
was also negative for segment 8 until day 2 pi when the 
LVI was detected in all 4 fish and HVI was detected in 
2/4 fish. Thereafter from day 3 pi, both viruses were de- 
tected in all organ samples. Statistical modelling indi- 
cates that the time course of the 2 strains, as indicated 
by detection of segment 8, differ (p < 0.05) for all 3 or- 
gans (Figure 1). This was most apparent in the early 
stages of the infection up to day 8 and especially notable 
in the gills (Figure 1A) during the first 4 days pi. Follow- 
ing day 8 pi, the LVI load almost stabilised in compari- 
son to the HVI load which continued to increase further 
up to day 15 pi. All uninfected control fish were negative 
for ISAV. 

RNA species specific analysis of segment 8 

The assays targeting RNA species corresponding specific- 
ally to the replication (cRNA) and transcription (mRNA) 
processes provided further evidence the LVI had success- 
fully entered and replicated in the gills early in infection 
(Table 1). The results also suggest HVI has a slower up- 
take and/or entry into the gills and therefore progresses 
more slowly into generalised infection. In the LVI group, 
mRNA was detected in all fish (gills) at day 1 pi and there- 
after until experimental end. In comparison, HVI mRNA 
was not first detected in gills until day 4 pi (in 3 fish) 
which correlated with the first large increase in HVI load 
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Figure 1 Statistical modelling of ISAV segment 8. Statistical modelling of viral segment 8 total RNA load profiles of high virulent (dashed line) 
and low virulent (solid line) strains measured by real-time RT-PCR from gill (A), anterior kidney (B) and heart (C) tissues post-immersion infection 
at Day 0. Dotted lines indicate 95% confidence intervals. 



by seg8 qPCR. Similarly, the LVI cRNA was also detected 
earlier in the gills than in the HVI group, although there 
was increased variation at the individual level in the latter. 
The difference in mRNA and cRNA production between 
strains was less pronounced in the kidney and the heart. 
In kidney, mRNA of both ISAV strains was detected at 
day 3 pi, however the number of positive individuals 
remained variable in the HVI group until day 7 pi. The 
first detection of cRNA in the kidney of the LVI group 
was in 1 fish on day 3 pi, followed by 3 fish on day 4 pi. 



The HVI cRNA was first detected in the kidney from 3 
fish on day 4 pi (on days 5 and 6 also). In heart, the first 
mRNA detections for both viruses were from single fish 
on days 2 and 4 respectively, prior to all 4 fish on days 4 
and 7 respectively. The first cRNA detections in heart oc- 
curred in 1 fish on day 3 pi in the LVI group and in 2 fish 
on day 4 in the HVI group. 

As expected, the Cp values indicated cRNA was present 
in lower amounts than mRNA [38,44]. Cross-checking Cp 
results with melt curve data indicated non-specificity 



McBeath et ah Veterinary Research 2014, 45:83 
http://www.veterinaryresearch.Org/content/45/1/83 



Page 5 of 13 



Table 1 RNA strand specific real-time RT-PCR 
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Number of fish (max. = 4) testing positive for mRNA or cRNA of low virulent (LVI) or highly virulent (HVI) virus using RNA strand specific real-time RT-PCR analysis 
of gill, anterior kidney and heart, na = not applicable. 



interfered with the mRNA assay when template levels 
were low (seg8 Cp>32), as discussed previously [38]. 
Contrary to this, the cRNA assay was more specific, even 
at very low levels (seg8 Cp > 38) (data not shown). 

Immune gene expression 

The expression of four immune markers, Type I and II 
IFN, Mx and yIP, was measured in gills, heart and anter- 
ior kidney to provide information on both innate and 
adaptive immune responses upon infection with the two 
different ISAV strains. All immune genes in infected fish 
were strongly up-regulated in comparison to the nega- 
tive control fish throughout the experiment (Table 2). 
The maximum increase in transcript levels for all four 
immune gene markers to both ISAV strains was higher 
in kidney and heart compared to gills. The highest max- 
imum increase was observed for Type II IFN in kidney. 
Statistical modelling indicates that the time course of ex- 
pression for all four genes differ between the two strains 
(/?<0.05) in gill (Figure 2), anterior kidney (Figure 3) 
and heart (Figure 4) indicating a differential immuno- 
logical host response to the two viruses. In all cases, the 
immune genes were stimulated significantly more in the 
LVI group than the HVI group up to day 7 or 8 pi. In 



contrast, after day 8, the immune genes were expressed 
at higher levels in the HVI fish. 

Segment 8 and immune gene correlation analysis 

The associations between the difference in virus segment 
8 of the low and the highly virulent viruses and the dif- 
ference in immune gene expression response to each 
strain were investigated using correlation. The rationale 
for this is that a positive association between an increas- 
ing difference in load between the two strains and an in- 
creasing immune response is indicative of a possible 
causal relationship between viral load and immune re- 
sponse. Association of the differences between the seg- 
ment 8 viral load and the expression of all immune 
genes in the gill and heart for both LVI and HVI 
(Table 3) were observed (p < 0.05), indicating an associ- 
ation between viral load as measured by segment 8 RNA 
quantity and immune gene expression in these organs. 

Segment 7 

Expression profiles for genomic segment 7 in gill, kidney 
and heart (Figure 5) were similar overall to those of seg- 
ment 8, although a higher quantity of segment 7 was 
produced by LVI than HVI up to day 8 pi. This was 



Table 2 Fold increase 


in expression 


of immune 


genes 












Gene 


Gill 






Kidney 






Heart 






LVI 


HVI 


Control 


LVI 


HVI 


Control 


LVI 


HVI 


Control 


Type 1 IFN 


7.7 


10.4 


2.7 


49.6 


24.6 


2.4 


23.6 


100.5 


1.8 


Mx 


71.6 


63 


1.4 


213.4 


100 


2.9 


257.2 


271.9 


2.6 


Type II IFN 


30.2 


47 


5.9 


1443.2 


485 


6.4 


394.5 


197.9 


7.6 


gIP 


93.2 


54 


2.3 


242.5 


255 


4.1 


290.7 


424.3 


9.8 



Fold increase from 6 h to the maximum peak of immune gene expression (on any given day) post-infection (pi) for fish infected with ISAV of either low (LVI) or 
high (HVI) virulence, or cell culture media to serve as negative controls. Figures in bold indicate the highest fold-induction for each gene in each organ. 
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Figure 2 Statistical modelling of immune genes in gill. Statistical modelling of 4 immune marker expression profiles, Type I IFN (A), Mx (B), 

Type II IFN (C) and ylP (D) in gills offish infected with high virulent (dashed line) or low virulent (solid line) strains measured by real-time RT-PCR. 
Dotted lines indicate 95% confidence intervals. 



more evident upon comparing a ratio of the two seg- 
ments for each strain (Additional file 2). Following day 8 
pi, unlike segment 8, there was little difference between 
the two viruses with regards to segment 7 RNA produc- 
tion as indicated by overlapping confidence intervals. In 
contrast to the results for segment 8, there was no evi- 
dence of an association in differences between the seg- 
ment 7 viral loads and differences in immune response 
for gill and heart, although positive correlation (p < 0.05) 
was observed in kidney (Table 3). 



Discussion 

In the present study, organs from fish infected with low 
virulent ISAV (LVI) exhibited a higher viral load in the 
first 8 days of the infection compared to the highly viru- 
lent virus (HVI). This was most prominent in the gills 
during the first 4 days post-infection. The rapid increase 
in LVI segment 8 RNA load suggests the occurrence of a 
fast primary replication phase in the gill. The RNA spe- 
cies specific assays provided further support for this, in- 
dicating the LVI successfully gained fast entry to host 
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Figure 3 Statistical modelling of immune genes in kidney. Statistical modelling of 4 immune marker expression profiles, Type I IFN (A), Mx 
(B), Type II IFN (C) and ylP (D) in anterior kidney of fish infected with high virulent (dashed line) or low virulent (solid line) strains measured by 
real-time PCR. Dotted lines indicate 95% confidence intervals. 



gill cells, and was both actively replicating and transcrib- 
ing in this early stage of infection. All four immune 
genes, representing both Type I and Type II IFN systems 
and two products thereof, were significantly differentially 
expressed when comparing the two ISAV isolates. This 
suggested differences in the host response towards the 
two virus infections. Statistical correlations provided 
further evidence suggesting a causal link between virus 
production (as measured by segment 8 RNA) and up- 



regulation of the studied immune genes in the gill and 
heart. 

In an immersion challenge, all mucosal surfaces, in- 
cluding gill, skin, eye and gut, are possible entry points 
for microbial infections. The gills have previously been 
suggested as the main entry port for ISAV [45]. However 
fin and skin have both been implicated as important 
entry points for viral haemorrhagic septicaemia (VHSV), 
infectious haematopoietic necrosis virus (IHNV) and koi 
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Figure 4 Statistical modelling of immune genes in heart. Statistical modelling of 4 immune marker expression profiles, Type I IFN (A), Mx (B), 

Type II IFN (C) and ylP (D) in heart offish infected with high virulent (dashed line) or low virulent (solid line) strains measured by real-time RT-PCR. 
Dotted lines indicate 95% confidence intervals. 



Table 3 Correlation analysis 



Segment 8 



Segment 7 



Immune gene 


Gill 


Kidney 


Heart 


Gill 


Kidney 


Heart 


Mx 


0.75** 


0.9*** 


0.9*** 


0.23 


0.72** 


0.12 


ylP 


0.61* 


0.43 


0.96*** 


-0.002 


0.46 


0.12 


Type I IFN 


0.61* 


0.43 


0.86*** 


-0.002 


0.56* 


-0.002 


Type II IFN 


0.82*** 


0.55* 


0.78** 


0.49 


0.56* 


0.2 



Correlation coefficients (r s ) for differences between strains with respect to virus 
burden and immune response. (*0.01 <p<0.05, **0.001 <p<0.01, ***p< 0.001). 



herpes virus (KHV) [46-48]. A cohabitation infection 
trial using the same HVI as the present study observed a 
relatively low level of infection in the gills, however few 
sampling points were used [34]. The currently presented 
viral load kinetics, backed up by the assays directly tar- 
geting cRNA and mRNA, suggest the presence of an 
early primary replication phase for the LVI in the gills, 
not observed for the HVI. Detection of both viruses in 
all fish at 6 h pi by seg8 real-time RT-PCR indicated 
rapid uptake of both viruses to the gills from the water. 
However, only the LVI amount displayed an early 
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Figure 5 Statistical modelling of ISAV segment 7. Statistical modelling of viral segment 7 total RNA load profiles of high virulent (dashed line) 
and low virulent (solid line) strains measured by real-time RT-PCR from gill (A), anterior kidney (B) and heart (C) tissues. Dotted lines indicate 95% 
confidence intervals. 



immediate increase which can initially be attributed to 
the rapid production of mRNA in the time points that 
follow. The replicative intermediate, cRNA, was also de- 
tected in some fish from day 1. As expected this was at 
much lower levels throughout the experiment [38,44]. 
Therefore the production of cRNA had a negligible ef- 
fect on the overall viral load compared to increasing 
mRNA production and accumulating genome (vRNA) 
from nascent virions. The subsequent effect of this early 
virus multiplication was clearly visible in the kidney and 



heart, where LVI was consistently detected in greater 
quantities each day up until day 8 pi. The early LVI rep- 
lication (as indicated by the specific detection of the 
replicative intermediate cRNA) and accumulation of 
nascent virus particles in the gills, prior to dissemin- 
ation, probably accounts for this. Interestingly, the puta- 
tively non-virulent ISAV HPRO is primarily associated 
with detection in the gills only [9,49]. The gill results 
suggest HVI progresses more slowly into a generalised 
infection. However, the fact that HVI may utilise an 
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alternative entry point not examined in this study, can- 
not be excluded. 

Virus virulence is a multifunctional trait and replica- 
tion efficiency is one factor that has been associated with 
increased virulence. A study into the differential viru- 
lence mechanisms of two strains of IHNV of high or low 
virulence in salmonid fish was shown to be linked with 
faster in vivo kinetics and replication of the highly viru- 
lent virus [50,51]. Early innate host immune response 
was demonstrated to play a critical role and it was pro- 
posed that the lag in immune system stimulation might 
give the faster replicating highly virulent virus enough 
time to reach a threshold and outrun the host antiviral 
response [51]. In addition to the replication advantage, it 
was also established that the highly virulent IHNV had 
an advantage at the entry stage of the infection cycle, 
when data from immersion and IP-injection experiments 
were compared [52]. Contrary to this, investigations on 
simian immunodeficiency virus (SIV) have suggested 
high replication rates of the natural virus do not corres- 
pond with increased virulence [53]. In the present study 
following immersion challenge, whilst the LVI rapidly 
gained entry to gill cells and replicated immediately, it 
ultimately reached a lower peak viral load than the HVI 
and caused a lower mortality. We could argue this pro- 
vides LVI with an evolutionary advantage over HVI. Rep- 
lication to relatively high levels without killing all the 
available hosts may allow the virus to persist in a popu- 
lation and potentially spread further. However, further 
information on virus shedding rates is required to inves- 
tigate this. Our data suggests that the LVI has an advan- 
tage at the entry stage in gills and also exhibited a 
higher replication rate during the early stages of infec- 
tion. This finding agrees with previous experimental data 
demonstrating slower replication of a high pathogenic 
ISAV compared to a low pathogenic ISAV in Atlantic 
Salmon Kidney (ASK) cells at 20 °C [35]. The reason for 
apparent faster replication of LVI at early time points re- 
mains unknown, but could be due to enhanced viral 
RNA-dependant RNA polymerase (RdRp) activity. This 
activity could also be in conjunction with differential 
properties of the viral glycoproteins, such as more effi- 
cient particle attachment or membrane fusion, as dem- 
onstrated for highly pathogenic influenza A viruses 
[54,55]. At present, functional information on ISAV 
RdRp's, and on the links between RdRp's and surface 
proteins is limited. 

The significant up-regulation of all four immune genes 
showed both LVI and HVI effectively stimulated the in- 
nate and adaptive immune systems. The response to 
both viruses was substantial and sustained. It was evident 
LVI disseminated more rapidly due to its replication at the 
earlier stages of infection, thus triggering an earlier sys- 
temic IFN-mediated immune response compared to that 



of HVI. Infection of TO cells with ISAV of low and high 
virulence suggested replication was a requirement for in- 
duction of immune genes [36]. The generation of dsRNA 
is a pivotal part of the replication process for many ssRNA 
viruses and longer stretches of dsRNA are key indicators 
of viral invasion to a cell [56]. In the present study, the de- 
tection of the replicative intermediary cRNA early in the 
gills signals its presence. Cells detect viral RNA and pro- 
teins via pathogen associated molecular pattern (PAMP) 
receptors, which in turn stimulate interferons and the 
antiviral response [32,57]. In this study, the more rapid 
systemic response induced by LVI might have provided a 
sufficient level of protection in a higher number of hosts, 
preventing this strain from reaching the damaging higher 
viral loads observed for HVI. In the latter half of the trial 
the immune genes of the HVI infected fish were up- 
regulated more in comparison to those of the LVI infected 
fish. Interferon responses to viral infection are usually 
transient and self-limited to avoid a prolonged anti-viral 
state which in itself can be detrimental to the host and 
interfere with haematopoiesis [57,58]. Indeed, the vast in- 
duction of cytokines and chemokines, generating a "cyto- 
kine storm" and overwhelming inflammatory responses, 
have been linked to highly pathogenic influenza virus 
pathogenesis [57,59,60]. The over- activation of IFN(3 and 
tumour necrosis factor-a (TNFa) creates a powerful pro- 
inflammatory response compared to that of low patho- 
genic influenza viruses, tipping the balance of the response 
towards inflammation, contributing to tissue damage [61]. 
The possibility that the increased mortality caused by 
HVI, which coincided with high expression of immune 
markers was caused by similar immune mechanisms 
should not be excluded. This study only focussed on a 
very small aspect of the immune response, therefore a 
more in-depth analysis of a greater number of immune re- 
sponse genes in immersion challenged fish would be ad- 
vantageous. What is clear from the present study, the 
immune response was sufficient to limit the infection by 
LVI while ineffectual at preventing HVI instigating a pro- 
gressive infection causing an eventual fatal outcome. 

RNA viruses have evolved diverse strategies to coun- 
teract and evade the host immune system. The influenza 
virus NS1 protein for example, is the primary antagonist 
of the innate immune response and remarkably effects 
many stages of the interferon response [62]. Two ISAV 
proteins, including the putative NS protein, have been 
linked to the antagonism of the IFN system [17,18], but 
it was not possible to relate these functions specifically 
here. In a study such as this, any viral IFN antagonistic 
effect within cells is likely to go unseen due to the vast 
cytokine inductions in neighbouring uninfected cells. 
Real-time PCR of segment 7 followed a similar profile to 
segment 8 in all three organs tested. However, correlation 
analysis suggested the segment 7 load, unlike segment 8, 



McBeath et ah Veterinary Research 2014, 45:83 
http://www.veterinaryresearch.Org/content/45/1/83 



Page 11 of 13 



is unlikely to be directly related to the immune gene ex- 
pression in either the gill or the heart. This has two pos- 
sible explanations. Firstly, the functions of the proteins 
encoded by these segments are very different and the ex- 
pression may be time dependent. In addition, the quan- 
tities of the required proteins are also substantially 
different and neither segment assay differentiates the two 
ORFs. The M protein is the most abundant protein in the 
virion [14], thus more is required compared to the non- 
or minor- structural proteins encoded by segment 7. 
Comparing the segment 7 to segment 8 ratios of both vi- 
ruses indicated the potentially interesting observation that 
HVI appears to generate less segment 7 RNA than LVI. 
The significance of this is unknown and remains an area 
of further investigation, although may partly explain the 
lack of correlation between segment 7 and immune genes. 

In conclusion, this study indicated that low virulent 
ISAV (LVI) replicated earlier in the gills and dissemi- 
nated throughout the host more efficiently following 
immersion challenge compared to the highly virulent 
virus (HVI). This suggests potential variation in tissue or 
even cell tropism. Rapid replication in gills stimulated a 
marked systemic immune response that may have pro- 
vided some protection against LVI-induced pathogen- 
esis, causing a limited infection. In contrast, the host 
response against HVI was less effectual, allowing the 
virus to reach a higher load causing a progressive infec- 
tion and ultimately inducing significantly more mortal- 
ity. Further work is required to elucidate the underlying 
factors that allow LVI to replicate more efficiently in the 
gills following immersion infection. Other factors such 
as variable HE activity allowing increased viral particle 
liberation to infect more cells, or differing abilities in im- 
mune system interference may have a role to play. In 
addition, because the HVI does not appear to access the 
gills as quickly, other entry points such as skin, eye, gut, 
should also be investigated. A greater understanding of 
the molecular and pathophysiological mechanisms of 
virulence and pathogenesis of LVI and HVI is critical in 
developing strategies and preventative measures to com- 
bat viral infection by ISAV. 

Additional files 



Competing interests 

The authors declare that they have no competing interests. 



Authors' contributions 

AM carried out and supervised all molecular experimental work, analysed 
results and drafted the manuscript YMH and AM carried out the RNA 
specific assays. MH and AM carried out statistical analysis. MH wrote the 
statistical sections in the manuscript. KF designed and organised the 
experimental challenge, prepared virus, and participated in the execution of 
the challenge and in data analysis. AM, MA, DHC, TM and IM participated in 
the design and execution of the experimental aquarium challenge and in 
data analysis. All authors contributed to the revision and editing of the 
manuscript and approved the final version. 

Acknowledgments 

The authors thank Mickael Fourrier and Rita Pettinello at Marine Scotland 
Science and Simon Weli at the Norwegian Veterinary Institute for laboratory 
assistance. Thanks also to the staff at the VESO Vikan research station, in 
particular our study director, Anne Berit Romstad. The research leading to 
these results has received funding from the European Community's Seventh 
Framework Programme (FP7, 2007-2013), Research Infrastructures action, 
under the grant agreement No. FP7-228394 (NADIR project), the Norwegian 
Research Council, HAVBRUK project No.186907 and 207024 and the Faroese 
research Council. 

Author details 

1 Marine Scotland Science, Marine Laboratory, Aberdeen, Scotland, UK. 
2 Norwegian Veterinary Institute, Oslo, Norway. 3 Food and Veterinary 
Authority, Torshavn, Faroe Islands. 

Received: 2 June 2014 Accepted: 24 July 2014 
Published: 21 August 2014 

References 

1. Aamelfot M, Dale OB, Falk K: Infectious salmon anaemia - pathogenesis 
and tropism. J Fish Dis 2014, 37:291-307. 

2. Thorud K, Djupvik HO: Infectious anaemia in Atlantic salmon (Salmo salar 
L). Bull Eur Assoc Fish Pathol 1 988, 8:1 09-1 1 1 . 

3. Bouchard D, Keleher W, Opitz HM, Blake S, Edwards KC, Nicholson BL: 
Isolation of infectious salmon anemia virus (ISAV) from Atlantic salmon 
in New Brunswick, Canada. Dis Aquat Organ 1999, 35:131-137. 

4. Lovely JE, Dannevig BH, Falk K, Hutchin L, MacKinnon AM, Melville KJ, 
Rimstad E, Griffiths SG: First identification of infectious salmon anaemia 
virus in North America with haemorrhagic kidney syndrome. Dis Aquat 
Organ 1999, 35:145-148. 

5. Mullins JE, Groman DB, Wadowska D: Infectious salmon anaemia in salt 
water Atlantic salmon (Salmo salar L.) in New Brunswick, Canada. Bull Eur 
Assoc Fish Pathol 1998, 18:1 10-1 14. 

6. Murray AG, Munro LA, Wallace IS, Berx B, Pendrey DJ, Fraser D, Raynard RS: 
Epidemiological investigation into the re-emergence and control of an 
outbreak of infectious salmon anaemia in the Shetland Islands, Scotland. 
Dis Aquat Organ 2010, 91:189-200. 

7. Rodger HD, Turnbull T, Muir F, Millar S, Richards RH: Infectious salmon 
anaemia (ISA) in the United Kingdom. Bull Eur Assoc Fish Pathol 1998, 
18:115-116. 

8. Anonymous: ISA hits the Faroes. Fish Farming Int 2000, 27:47. 

9. Christiansen DH, 0stergaard PS, Snow M, Dale OB, Falk K: A low-pathogenic 
variant of infectious salmon anemia virus (ISAV-HPR0) is highly prevalent 
and causes a non-clinical transient infection in farmed Atlantic salmon 
(Salmo salar L.) in the Faroe Islands. J Gen Virol 201 1, 92:909-918. 

10. Godoy MG, Aedo A, Kibenge MJ, Groman DB, Yason CV, Grothusen H, 
Lisperguer A, Calbucura M, Avendano F, Imilan M, Jarpa M, Kibenge FS: First 
detection, isolation and molecular characterization of infectious salmon 
anaemia virus associated with clinical disease in farmed Atlantic salmon 
(Salmo salar) in Chile. BMC Vet Res 2008, 4:28. 

1 1 . Rimstad E, Dale OB, Dannevig BH, Falk K: Infectious Salmon Anaemia. In Fish 
Diseases and Disorders, Volume 3: Viral, Bacterial and Fungal Infections. Edited by 
Woo PTK, Bruno DW. Oxfordshire: CAB International; 201 1:143-165. 

12. Mjaaland S, Rimstad E, Falk K, Dannevig BH: Genomic characterization of 
the virus causing infectious salmon anemia in Atlantic salmon (Salmo 
salar L): an orthomyxo-like virus in a teleost. J Virol 1997, 71:7681-7686. 

13. Aspehaug V, Mikalsen AB, Snow M, Biering E, Villoing S: Characterization of 
the infectious salmon anemia virus fusion protein. J Virol 2005, 
79:12544-12553. 



Additional file 1: Primers and probes used in RT-qPCR and RNA 
specific RT-qPCR. Table of details relating to all primer and probe 
sequences used in the standard real-time RT-PCR analysis and the RNA 
specific real-time RT-PCR analysis. 

Additional file 2: Difference in expression of ISAV segment 7 
compared to segment 8. Expression ratio of segment 7 compared to 
segment 8 in gill (A), anterior kidney (B) and heart (C) in fish infected 
with either low virulent (LVI) or highly virulent (HVI) ISAV. 



McBeath et ah Veterinary Research 2014, 45:83 
http://www.veterinaryresearch.Org/content/45/1/83 



14. Falk K, Aspehaug V, Vlasak R, Endresen C: Identification and 
characterization of viral structural proteins of infectious salmon anemia 
virus. J Virol 2004, 78:3063-3071. 

15. Rimstad E, Mjaaland S, Snow M, Mikalsen AB, Cunningham CO: 
Characterization of the infectious salmon anemia virus genomic 
segment that encodes the putative hemagglutinin. J Virol 2001, 
75:5352-5356. 

16. Biering E, Falk K, Hoel E, Thevarajan J, Joerink M, Nylund A, Endresen C, 
Krossoy B: Segment 8 encodes a structural protein of infectious salmon 
anaemia virus (ISAV); the co-linear transcript from Segment 7 probably 
encodes a non-structural or minor structural protein. Dis Aquat Organ 
2002, 49:117-122. 

1 7. Garcia-Rosado E, Markussen T, Kileng O, Baekkevold ES, Robertsen B, 
Mjaaland S, Rimstad E: Molecular and functional characterization of two 
infectious salmon anaemia virus (ISAV) proteins with type I interferon 
antagonizing activity. Virus Res 2008, 133:228-238. 

18. McBeath AJA, Collet B, Paley R, Duraffour S, Aspehaug V, Biering E, 
Secombes CJ, Snow M: Identification of an interferon antagonist protein 
encoded by segment 7 of infectious salmon anaemia virus. Virus Res 
2006, 115:176-184. 

19. Ramly RB, Olsen CM, Braaen S, Rimstad E: Infectious salmon anaemia virus 
nuclear export protein is encoded by a spliced gene product of 
genomic segment 7. Virus Res 2013, 177:1-10. 

20. Falk K, Namork E, Rimstad E, Mjaaland S, Dannevig BH: Characterization of 
infectious salmon anemia virus, an orthomyxo-like virus isolated from 
Atlantic salmon (Saimo salar L). J Virol 1997, 71:9016-9023. 

21. Snow M, McKay P, McBeath AJ, Black J, Doig F, Kerr R, Cunningham CO, 
Nylund A, Devoid M: Development, application and validation of a 
Taqman real-time RT-PCR assay for the detection of infectious salmon 
anaemia virus (ISAV) in Atlantic salmon (Saimo salar). Dev Biol (Basel) 
2006, 126:133-145. 

22. Dannevig B, Falk K, Skjerve E: Infectivity of internal tissues of Atlantic 
salmon, Saimo salar L, experimentally infected with the etiologic agent 
of infectious salmon anemia (ISA). J Fish Dis 1994, 17:613-622. 

23. Hetland DL, Dale OB, Skjodt K, Press CM, Falk K: Depletion of CD8 alpha 
cells from tissues of Atlantic salmon during the early stages of infection 
with high or low virulent strains of infectious salmon anaemia virus 
(ISAV). Dev Cornp Immunol 2011, 35:817-826. 

24. Kibenge FS, Kibenge MJ, Groman D, McGeachy S: In vivo correlates of 
infectious salmon anemia virus pathogenesis in fish. J Gen Virol 2006, 
87:2645-2652. 

25. Mjaaland S, Markussen T, Sindre H, Kj0glum S, Dannevig BH, Larsen S, 
Grimholt U: Susceptibility and immune responses following experimental 
infection of MHC compatible Atlantic salmon (Saimo salar L.) with 
different infectious salmon anaemia virus isolates. Arch Virol 2005, 
150:2195-2216. 

26. Ritchie RJ, McDonald JT, Glebe B, Young-Lai W, Johnsen E, Gagne N: 
Comparative virulence of Infectious salmon anaemia virus isolates in 
Atlantic salmon, Saimo salar L. J Fish Dis 2009, 32:1 57-171. 

27. Markussen T, Jonassen CM, Numanovic S, Braaen S, Hjortaas M, Nilsen H, 
Mjaaland S: Evolutionary mechanisms involved in the virulence of 
infectious salmon anaemia virus (ISAV), a piscine orthomyxovirus. 
Virology 2008, 374:515-527. 

28. Mjaaland S, Hungnes O, Teig A, Dannevig BH, Thorud K, Rimstad E: 
Polymorphism in the infectious salmon anemia virus hemagglutinin 
gene: importance and possible implications for evolution and ecology of 
infectious salmon anemia disease. Virology 2002, 304:379-391. 

29. Markussen T, Sindre H, Jonassen CM, Tengs T, Kristoffersen AB, Ramsell J, 
Numanovic S, Hjortaas MJ, Christiansen DH, Dale OB, Falk K: Ultra-deep 
pyrosequencing of partial surface protein genes from infectious Salmon 
Anaemia Virus (ISAV) suggest novel mechanisms involved in transition 
to virulence. PLoS One 2013, 8:e81 571 . 

30. Langevin C, Aleksejeva E, Passoni G, Palha N, Levraud JP, Boudinot P: The 
antiviral innate immune response in fish: evolution and conservation of 
the IFN system. J Mol Biol 201 3, 425:4904-4920. 

31. Zou J, Secombes CJ: Teleost fish interferons and their role in immunity. 
DevComp Immunol 2011, 35:1376-1387. 

32. Collet B: Innate immune responses of salmonid fish to viral infections. 
Dev Comp Immunol 2014, 43:160-173. 

33. McBeath AJA, Snow M, Secombes CJ, Ellis AE, Collet B: Expression kinetics 
of interferon and interferon-induced genes in Atlantic salmon (Saimo 



Page 12 of 13 



salar) following infection with infectious pancreatic necrosis virus and 
infectious salmon anaemia virus. Fish Shellfish Immunol 2007, 22:230-241. 

34. Jorgensen SM, Hetland DL, Press CM, Grimholt U, Gjoen T: Effect of early 
infectious salmon anaemia virus (ISAV) infection on expression of MHC 
pathway genes and type I and II interferon in Atlantic salmon (Saimo 
salar L.) tissues. Fish Shellfish Immunol 2007, 23:576-588. 

35. Svingerud T, Holand JK, Robertsen B: Infectious salmon anemia virus 
(ISAV) replication is transiently inhibited by Atlantic salmon type I 
interferon in cell culture. Virus Res 2013, 177:163-170. 

36. Workenhe ST, Hori TS, Rise ML, Kibenge MJ, Kibenge FS: Infectious salmon 
anaemia virus (ISAV) isolates induce distinct gene expression responses 
in the Atlantic salmon (Saimo salar) macrophage/dendritic-like cell line 
TO, assessed using genomic techniques. Mol Immunol 2009, 
46:2955-2974. 

37. McBeath A, Aamelfot M, Christiansen DH, Matejusova I, Markussen T, 
Kaldhusdal M, Dale OB, Weli SC, Falk K: Immersion challenge with low and 
highly virulent infectious salmon anaemia virus reveals different 
pathogenesis in Atlantic salmon, Saimo salar L. J Fish Dis. in press. 

38. McBeath A, Bain N, Fourrier M, Collet B, Snow M: A strand specific real- 
time RT-PCR method for the targeted detection of the three species 
(vRNA, cRNA and mRNA) of infectious salmon anaemia virus (ISAV) 
replicative RNA. J Virol Methods 2013, 187:65-71. 

39. Zou J, Carrington A, Collet B, Dijkstra JM, Yoshiura Y, Bols N, Secombes C: 
Identification and bioactivities of IFN-gamma in rainbow trout 
Oncorhynchus mykiss: the first Thl-type cytokine characterized 
functionally in fish. J Immunol 2005, 175:2484-2494. 

40. R Core Team: ft A language and environment for statistical computing. 
Vienna, Austria: R Foundation for Statistical Computing; [http://www.R- 
project.org/]. (Accessed 29 November 2012). 

41. Cleveland WS: Robust locally weighted regression and smoothing 
scatterpoints. J Am Stat Assoc 1979, 74:829-836. 

42. Searle SR: Linear Models. New York: John Wiley and Sons; 1 971 . 

43. Spearman C: The proof and measurement of association between two 
things. Am J Psychol 1904, 15:72-101. 

44. Hatada E, Hasegawa M, Mukaigawa J, Shimizu K, Fukuda R: Control of 
influenza virus gene expression: quantitative analysis of each viral RNA 
species in infected cells. J Biochem 1989, 105:537-546. 

45. Mikalsen AB, Teig A, Helleman AL, Mjaaland S, Rimstad E: Detection of 
infectious salmon anaemia virus (ISAV) by RT-PCR after cohabitant 
exposure in Atlantic salmon Saimo salar. Dis Aquat Organ 2001, 47:175-181. 

46. Costes B, Raj VS, Michel B, Fournier G, Thirion M, Gillet L, Mast J, Lieffrig F, 
Bremont M, Vanderplasschen A: The major portal of entry of koi 
herpesvirus in Cyprinus carpio is the skin. J Virol 2009, 83:2819-2830. 

47. Harmache A, LeBerre M, Droineau S, Giovannini M, Bremont M: 
Bioluminescence imaging of live infected salmonids reveals that the fin 
bases are the major portal of entry for Novirhabdovirus. J Virol 2006, 
80:3655-3659. 

48. Quillet E, Dorson M, Aubard G, Torhy C: In vitro assay to select rainbow 
trout with variable resistance/susceptibility to viral haemorrhagic 
septicaemia virus. Dis Aquat Organ 2007, 76:7-16. 

49. McBeath AJ, Bain N, Snow M: Surveillance for infectious salmon anaemia 
virus HPR0 in marine Atlantic salmon farms across Scotland. Dis Aquat 
Organ 2009, 87:161-169. 

50. Penaranda MM, Purcell MK, Kurath G: Differential virulence mechanisms of 
infectious hematopoietic necrosis virus in rainbow trout (Oncorhynchus 
mykiss) include host entry and virus replication kinetics. J Gen Virol 2009, 
90:2172-2182. 

51 . Purcell MK, Garver KA, Conway C, Elliott DG, Kurath G: Infectious 
haematopoietic necrosis virus genogroup-specific virulence mechanisms 
in sockeye salmon, Oncorhynchus nerka (Walbaum), from Redfish Lake, 
Idaho. J Fish Dis 2009, 32:61 9-631 . 

52. Wargo AR, Kurath G: In vivo fitness associated with high virulence in a 
vertebrate virus is a complex trait regulated by host entry, replication, 
and shedding. J Virol 201 1, 85:3959-3967. 

53. Pandrea I, Sodora DL, Silvestri G, Apetrei C: Into the wild: simian 
immunodeficiency virus (SIV) infection in natural hosts. Trends Immunol 
2008, 29:419-428. 

54. Pappas C, Aguilar PV, Basler CF, Solorzano A, Zeng H, Perrone LA, Palese P, 
Garcia-Sastre A, Katz JM, Tumpey TM: Single gene reassortants identify a 
critical role for PB1, HA, and NA in the high virulence of the 1918 
pandemic influenza virus. Proc Natl Acad Sci USA 2008, 1 05:3064-3069. 



McBeath et al. Veterinary Research 2014, 45:83 
http://www.veterinaryresearch.Org/content/45/1/83 



Page 13 of 13 



55. Grimm D, Staeheli P, Hufbauer M, Koerner I, Martinez-Sobrido L, Solorzano 
A, Garcia-Sastre A, Haller 0, Kochs G: Replication fitness determines high 
virulence of influenza A virus in mice carrying functional Mxl resistance 
gene. Proc Natl Acad Sci USA 2007, 1 04:6806-681 1 . 

56. Wang Q, Carmichael GG: Effects of length and location on the cellular 
response to double-stranded RNA. Microbiol Mol Biol Rev 2004, 
68:432-452. 

57. Kash JC: Applications of high-throughput genomics to antiviral research: 
evasion of antiviral responses and activation of inflammation during 
fulminant RNA virus infection. Antiv Res 2009, 83:10-20. 

58. Lin Q, Dong C, Cooper MD: Impairment of T and B cell development by 
treatment with a type I interferon. J Exp Med 1 998, 1 87:79-87. 

59. Cilloniz C, Shinya K, Peng X, Korth MJ, Proll SC, Aicher LD, Carter VS, Chang 
JH, Kobasa D, Feldmann F, Strong JE, Feldmann H, Kawaoka Y, Katze MG: 
Lethal influenza virus infection in macaques is associated with early 
dysregulation of inflammatory related genes. PLoS Pathog 2009, 

5:e 1000604. 

60. Marsolais D, Hahm B, Walsh KB, Edelmann KH, McGavern D, Hatta Y, 
Kawaoka Y, Rosen H, Oldstone MB: A critical role for the sphingosine 
analog AAL-R in dampening the cytokine response during influenza 
virus infection. Proc Natl Acad Sci USA 2009, 1 06:1 560-1 565. 

61. Lee SM, Gardy JL, Cheung CY, Cheung TK, Hui KP, Ip NY, Guan Y, Hancock 
RE, Peiris JS: Systems-level comparison of host-responses elicited by avian 
H5N1 and seasonal H1N1 influenza viruses in primary human 
macrophages. PLoS One 2009, 4:e8072. 

62. Wang BX, Fish EN: The yin and yang of viruses and interferons. Trends 
Immunol 2012, 33:190-197. 



doi:10.1 186/S13567-014-0083-X 

Cite this article as: McBeath ef al:. Low virulent infectious salmon 
anaemia virus (ISAV) replicates and initiates the immune response 
earlier than a highly virulent virus in Atlantic salmon gills. Veterinary 
Research 2014 45:83. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at /^'\ n! _, ul _-| r Q r,tr=l 

www.biomedcentral.com/submit ammBa central 



V 



